Diverging probability density functions for flat-top solitary waves 
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We investigate the statistics of flat-top solitary wave parameters in the presence of weak multi- 
plicative dissipative disorder. We consider first propagation of solitary waves of the cubic-quintic 
nonlinear Schrodinger equation (CQNLSE) in the presence of disorder in the cubic nonlinear gain. 
We show by a perturbative analytic calculation and by Monte Carlo simulations that the probabil- 
ity density function (PDF) of the amplitude r\ exhibits loglognormal divergence near the maximum 
possible amplitude r] m , a behavior that is similar to the one observed earlier for disorder in the linear 
gain [A. Peleg et al., Phys. Rev. E 72, 027203 (2005)]. We relate the loglognormal divergence of 
the amplitude PDF to the super-exponential approach of n to rj m in the corresponding deterministic 
model with linear/nonlinear gain. Furthermore, for solitary waves of the derivative CQNLSE with 
weak disorder in the linear gain both the amplitude and the group velocity /3 become random. We 
therefore study analytically and by Monte Carlo simulations the PDF of the parameter p, where 
p = T)/{1 — e s [3/2) and e s is the self-steepening coefficient. Our analytic calculations and numerical 
simulations show that the PDF of p is loglognormally divergent near the maximum p- value. 

PACS numbers: 05.45.Yv, 05.40.-a, 47.54.-r, 42.65.Tg 
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I. INTRODUCTION 

Flat-top solitary waves are coherent patterns, which 
exist as a result of a balance between disper- 
sion/diffraction and competing nonlinearities, where the 
low order nonlinearity is "focusing" while the high order 
nonlinearity is "defocusing" [l|, 0, IH . When the intensity 
of the field is relatively small, the low order nonlinear- 
ity is dominant, and consequently, the solitary waves are 
narrow and have a shape that is similar to that of conven- 
tional solitons. However, when the intensity increases, 
the high order nonlinearity becomes dominant and leads 
to the broadening of the pulse shape and to the gen- 
eration of a typical table-top pattern. Flat-top solitary 
waves appear as solutions to nonlinear wave models in 
many areas of physics, including nonlinear optics yj, y, [5[ 
fluid dynamics d, 0], and plasma physics 0, 0- As a 
result, they have been the subject of intensive research 
efforts in recent years. The interest in flat-top solitary 
waves is further enhanced since they are used in pat- 
tern formation theory to explain the emergence of fronts 
(kinks) from localized coherent structures such as soli- 
tons and solitary waves [E EE] ■ Many of the systems in 
which flat-top solitary waves appear can be influenced 
by processes involving noise or disorder. When the dis- 
order is strong the solitary waves are usually destroyed, 
whereas, when it is weak, the solitary waves can form 
and evolve. In the latter case one is mainly concerned 
with the statistics of the solitary wave parameters. 

In this study we focus attention on an important type 
of disorder, which we call multiplicative dissipative dis- 
order. This type of disorder is characterized by the fol- 
lowing properties: (1) the disorder affects the amplitude 
of the solitary wave in first order of the perturbation; (2) 
the disorder term in the nonlinear wave model is mul- 



tiplicative. Dissipative disorder can appear in systems 
described by nonlinear wave equations in various forms. 
Two of the most common forms are disorder in the linear 
gain/loss coefficient and disorder in the cubic nonlinear 
gain/loss coefficient. Disorder in the linear gain coeffi- 
cient can appear in optical fiber communication systems 
due to randomness in the gain of amplifiers that are posi- 
tioned along the fiber line to compensate for the loss [ll| . 
Moreover, such disorder appears in massive multichannel 
optical fiber transmission as a result of the interplay be- 
tween Raman-induced energy exchan ge i n pulse collisions 
and bit pattern randomness EE EE HI], EE GE] • ^ n this 
case, the disorder can lead to relatively high bit-error- 
rate values and intermittent dynamics of pulse param- 
eters [IE GE GH- We point out that in all the studies 
reported in Refs. [Il|, GJ, GE GI GE EE EE El weak 
disorder was considered. In addition, both linear and 
cubic nonlinear disorder in the gain can emerge in an 
active nonlinear medium due to random variations with 
distance in the linear/nonlinear gain/loss coefficient. 

We consider two nonlinear wave models, which possess 
flat-top solitary wave solutions: the cubic-quintic non- 
linear Schrodinger equation (CQNLSE) and the deriva- 
tive CQNLSE (DCQNLSE). A third model, the ex- 
tended Korteweg-de Vries (eKdV) equation, is also 
briefly discussed. The CQNLSE is a simple nonintc- 
grable extension of the cubic nonlinear Schrodinger equa- 
tion (CNLSE) possessing solitary wave solutions. The 
CQNLSE describes a variety of physical systems includ- 
ing pulse propa gati on in semiconductor-doped optical 
fibers ELJj, llS. |20L l2ll [22l . |23| . laser-plasma interac- 
tion [7u24l|. and Bose-Einstein condensates [IE HE HE 
[28j . Another important reason for the interest in the 
CQNLSE is that due to its nonintegrability it allows one 
to observe dynamical effects that do not exist in the 
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CNLSE, e.g., emission of continuous radiation in two- 
soliton collisions [2!|[3(|. Furthermore, the cubic-quintic 
complex Ginzburg-Landau equation, which is a gener- 
alization of the CQNLSE including dissipative terms, is 
known to describe even a wider range of physical systems, 
including, for ex ample, convection and pattern formation 
in fluids [IJlEM H M 13 and mode-locked lasers 

SlIlSTHlii 

The derivative CNLSE is an extension of the CNLSE 
that, in the context of nonlinear optics, takes into ac- 
count the effects of self-steepening [19|, |39|, H(| E| ■ For 
short optical pulses propagating in semiconductor-doped 
fibers both quintic nonlinearity and self-steepening are 
important. Consequently, one expects that the derivative 
CQNLSE, which takes into account both effects, would 
provide a more accurate description of the propagation in 
this case. We note that a variant of the derivative CNLSE 
is known to describe propaga tion of Alfven waves in mag- 
netized plasmas [HI, l43l. 144]. Moreover, it was recently 
shown theoretically and experimentally that a variant of 
the DCQNLSE accurately describes propagation of high- 
intensity pulses in cascaded-quadratic nonlinear media 

The eKdV equation, which is also known as the Gard- 
ner equation, is an intcgrable model that describes inter- 
facial waves in a two-layer system |46l l47l |48| as well as 
stratified shear flow in the ocean 4&5oj]. It provides 
a possible explanation to observations of large- amp litude 
fiat-top solitary waves in coastal zones 
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In a previous work (54[ we studied the effects of weak 
disorder in the linear gain coefficient on solitary waves 
of the CQNLSE. We showed analytically (by employing 
an adiabatic perturbation method) and by Monte Carlo 
simulations that the probability density function (PDF) 
of the solitary wave amplitude has a loglognormal diverg- 
ing form in the vicinity of the maximum possible ampli- 
tude. Since solitary waves with amplitude values close to 
the maximum possible amplitude have a table-top shape, 
this finding means that the amplitude PDF of flat-top 
solitary waves exhibits loglognormal divergence. We also 
conjectured that similar loglognormal divergence should 
be observed for disorder in the nonlinear gain. However, 
the full analytic form of the amplitude PDF for disor- 
der in the nonlinear gain was not obtained and the con- 
jecture was not tested by numerical simulations. Thus, 
the important question concerning the generality of the 
loglognormal divergence of the amplitude PDF for flat- 
top solitary waves was not fully answered, even for the 
CQNLSE. Furthermore, it is not clear whether loglog- 
normal divergence can be observed in other nonlinear 
wave models possessing solitary wave solutions. In the 
current paper we address these questions in detail. Wc 
start by considering propagation of solitary waves of the 
CQNLSE in the presence of weak disorder in the cubic 
nonlinear gain. The case of disorder in the cubic nonlin- 
ear gain is particularly important since cubic gain/loss 
is very common in many systems described by the com- 
plex Ginzburg-Landau equation (see, for example, Ref. 



[HI and references therein). We calculate the amplitude 
PDF analytically by employing an adiabatic perturba- 
tion method and validate its loglognormal divergence by 
Monte Carlo simulations. We then turn to study propa- 
gation of solitary waves of the DCQNLSE in the presence 
of weak disorder in the linear gain. In this case both the 
amplitude and the group velocity randomly vary during 
propagation. We therefore study the PDF of a new pa- 
rameter, which is the ratio between the amplitude and 
a linear function of the group velocity. We find that the 
PDF of this new parameter is loglognormally divergent 
near the parameter's maximum value. We conclude by 
a brief discussion of the dynamic mechanism responsible 
for the loglognormal divergence of the PDFs, and of the 
possibility to observe similar statistical behavior in the 
eKdV model. 

The material in the rest of the paper is organized as fol- 
lows. In Sec. UH we study propagation of solitary waves 
of the CQNLSE in the presence of disorder in the cubic 
nonlinear gain. In Sec. IIIH we investigate propagation of 
solitary waves of the DCQNLSE in the presence of dis- 
order in the linear gain. Sections IIVI and [V] are reserved 
for discussion and conclusions, respectively. In Appendix 
|A1 we describe a method for identifying loglognormal di- 
vergence in numerical data. Finally, Appendixes [B] and 
[Cl are devoted to calculation of the PDF by employing 
the Fokker-Planck approach and by working within Ito's 
interpretation. 



II. CUBIC-QUINTIC NLSE WITH DISORDER 
IN THE CUBIC GAIN/LOSS COEFFICIENT 

Consider the dynamics described by the CQNLSE with 
disorder in the cubic gain/loss coefficient: 

id z i> + 8%tl> + 2|VfV - e,MV = ie£(z)\iP\ 2 TP, (1) 

where the disorder £(z) is zero in average and short cor- 
related in z: 



<£(*)> = 0, D8(z - z'). 



(2) 



In the context of nonlinear optical waveguides ip is pro- 
portional to the envelope of the electric field, z is the 
propagation distance, t is a retarded time, e q is the quin- 
tic nonlinearity coefficient, < e <C 1 is the cubic gain 
coefficient and D is the disorder intensity. The terms 
e 9 |V'| 4 V' an d * e £(- 2 OI'0| 2 V' describe the effects of quintic 
nonlinearity and disorder in the cubic nonlinear gain/loss 
coefficient, respectively. When e = 0, Eq. (fT]) sup- 
ports stable solitary wave solutions of the form [551 ]: 
ip„(t, z) = * s (x) exp(ix), where 



*.(a0 



[(1 - V 2 /Vm) 1/2 cosh(2x) + l 



1/2 ' 



(3) 



r] m = (4e,/3)- 1 / 2 , x = a + P{t-y) + {rf - f3 2 )z, and 
x = rj(t — y — 2f3z). In these relations the parameters 
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r/, p, a, y are related to the amplitude, frequency, phase 
and position of the solitary wave, respectively. Note that 
the solitary wave solution ip a exists provided that i] < rj m . 

We study the dynamics of the solitary wave ips as de- 
scribed by Eq. (fT]). Since we are interested in flat-top 
solitary waves we focus attention on the case e q > 0. We 
also assume that 4De 2 z < 1, so that for most of the 
disorder realizations the dynamics of the solitary wave 
amplitude is not yet influenced by the 0(e 2 ) radiation 
instability effects [9, Hi| . The dynamics of the parameter 
77 is obtained by using energy balance considerations: 



dt\^\ 2 = 2e£(z) 



dt\i/>\ 4 



(4) 




In order to solve Eq. ^ we employ the adiabatic per- 
turbation method, which has been extensively used in 
revious studies of the CQNLSE, see, for example, Ref. 
This calculation yields: 



d_ 

dz 



arctanh ( 

7?m 



arctanh ( ) — 



Vn 



Vn 



(5) 



Furthermore, within the framework of Stratonovich's in- 
terpretation [53, \M, um E q- © 

can be transformed into 



drj 

dz 



= 4e£(z) (vl - V 2 ) 



r) m arctanh 



V 



Changing variables from 77 to v, where 

dr) 



I 



(6) 



(7) 



{Vm ~ V 2 ) Vim arctanh (r//r? m ) - rj\ ' 
we obtain the equation dv/dz = 4e£(z), whose solution 



FIG. 1: The probability density function F c {rj) at z = 10 for 
D = 3, e q = 0.5, e = 0.03 and t?(0) = 1. The solid curve 
corresponds to the analytic result obtained by using Eqs. (jSJ) 
and (|10|) . The squares represent the result of Monte Carlo 
simulations with Eq. (TTJ , while the circles stand for the result 
of Monte Carlo simulations with Eq. |(SJ . 
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v(z)-v(0) = 
= AeX(z), 



drj 



t?(0) (V m - V 2 ) [??marctanh (r)/r]m) - rj\ 



(8) 



where X(z) — L dz'£(z'). Notice that X(z) can be re- 
garded as a sum over many independent random vari- 
ables. Consequently, according to the central limit theo- 
rem its PDF approaches a Gaussian PDF of the form 



F(X) = (27rL>z)" 1/2 cxp [-X 2 /(2Dz)] 



(9) 



Equation |8]) defines a monotonously increasing function 
X = q(j]) on < r\ < r\ m . Changing variable from X to 
T] while employing Eqs. © and §§§ we obtain that the 
PDF of rj is given by 



(327rL>e 2 z)- 1 /2 exp [-f( v )/(2Dz)] 
{"nli " V 2 ) [Vm arctanh (77/77™) - 77] 



(10) 



for < 77 < ?7 m and F c (rj) = elsewhere. To calculate 
the value of F c (rj) for a given 77 one numerically solves Eq. 



FIG. 2: Blowup of the data in Fig. [T] in the vicinity of 77„ 
showing the divergence of F c {rj), 



© for X = q{rj) and then substitutes the result into Eq. 
([TO")) . The graph of F c (rj) obtained by this calculation is 
shown in Figs. Q]and[21 

As can be seen from Figs. [T]and[2l F c (rf) diverges in 
the vicinity ofr] m . In order to characterize the divergence 
we obtain an approximate analytic solution of Eq. (jHJ) for 
77 near ?7 m . We first write the integral in Eq. (jHJ) as a sum 
of two integrals: 



4eX(z) = \ 

J V (0) 



drj 



(o) Wm - V 2 ) [Vm arctanh {rj/r] m ) - rj] 

dv 

v (z) {"nli - V 2 ) [?7marctanh (77/77™) - 77] ' 



11) 



where r)(z) is a constant satisfying 77(2) < rj(z) < r) mi and 
both 77(5) and rj(z) are close to ?7 m . Since both limits 
of the first integral on the right hand side of Eq. (fTTT) 
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are constants this integral is a constant that we denote 
by c\. We also denote br\ = rj m — r\ and notice that 
< Sr](z)/r] m < Sr](z)/r] m <C 1. We can therefore expand 
the integrand in the second integral on the right hand side 
of Eq. (jlip about Srj = 0, keeping terms up to order 5r) 
in the denominator. This calculation yields 



AeX(z) ~ ci + 4" / 



Sr](z) 



dSl] 



(12) 



Integrating over Sr/ we arrive at 

' e 2 5rj(z) 



X(z) ~ In 



hi 



2?7„ 



/cW(4e^), (13) 



where c is another constant. Since the normally dis- 
tributed random variable X(z) is related to Srj(z) via 
a double logarithm, we say that 8rj{z) is loglognormally 
distributed. Using Eqs. (|SJ) and (fl"3|) . and changing vari- 
ables from X to Srj we obtain 



In 



e 2 5r) 

2?7m 



xexp ■ 



ln 2 [-ln[(e 2 ^)/(2r/ m )] /c] 

32^^ 



(14) 



from which it follows that the divergence of F c (rf) near 
?7 m is loglognormal. 

To check the analytic predictions given by Eqs. (JTUJ) 
and (114[) we performed Monte Carlo simulations with Eq. 
(TIJ with 1.09 x 10 5 disorder realizations. We considered 
the parameter values D = 3, e q = 0.5 (corresponding to 
r) m ~ 1.22474) and e = 0.03. The initial condition was 
taken in the form of the solitary wave ip s with 77(0) = 1, 
/3(0) = 0, y(0) = 0, and a(0) = 0. We carried out the 
simulations up to a final distance = 10, for which 
the disorder strength is 4De 2 Zf = 0.3. Equation (TIJ was 
integrated by employing a split-step method that is of 
fourth order with respect to the z-step dz (60j . The lin- 
ear part id z ip = —d 2 ip was advanced efficiently via an 
evaluation of the operator exponential in Fourier space 
and the nonlinear part id z ip — Sq^^ip+^e^z) — 2)|?/>| 2 ?/> 
was advanced via a fourth order Runge-Kutta scheme. 
To overcome numerical errors resulting from radiation 
emission and the use of periodic boundary conditions we 
applied the method of artificial damping in the vicinity 
of the boundaries of the computational domain. (See 
Refs. [13, HH [6l[ for other examples where the same 
method was successfully used). The size of the domain 
was taken to be — L < t < L with L = 16ir so that the ab- 
sorbing layers do not affect the dynamics of the solitary 
waves. The i-step and z-step were taken as At = 0.01 
and Az = 0.001, respectively. 

Figure [T] shows the 77-PDF obtained by the simulations 
as well as the analytic prediction obtained with Eqs. © 
and (fiTj| and the PDF obtained by Monte Carlo simula- 
tions with Eq. (|S]). Figure shows a blowup of the same 




FIG. 3: Gc(Sri) vs gc(Sri) for the same parameters considered 
in Figs. [T]and[2] The triangles represent the analytic predic- 
tion of Eqs. © and (|l(jp . the squares stand for the result of 
numerical simulations with Eq. {TJ and the circles correspond 
to the result obtained by Monte Carlo simulations with Eq. 
p]). The solid line is a linear fit of the squares with a slope 
of 0.97. 



data in the neighborhood of rj m . The good agreement 
between the three results strongly indicates that the di- 
vergence is indeed loglognormal. In Fig. Owe present a 
more sensitive analysis of this divergence that is based 
on the procedure described in Appendix [S] for detect- 
ing loglognormal divergence in numerical data. Follow- 
ing this procedure we plot G c (5t]) versus g c ($v)i where 
G c (Sr]) and g c (Sr]) are defined by Eqs. (|A"4|l and (|A5|) . 
respectively. It is seen that the data obtained by nu- 
merical simulations with Eq. (fT]) lies on a straight line 
with a slope 0.97, which is very close to the theoretically 
predicted value of 1. Therefore, this analysis provides 
further support in favor of the loglognormal divergence 
of F c {rj). Combining this observation with the result of 
Ref. [541], we conclude that both disorder in the linear 
gain/loss coefficient and disorder in the cubic nonlinear 
gain/loss coefficient lead to the same type of divergence 
of the 77-PDF. 



III. MULTIPLICATIVE-DISSIPATIVE 
DISORDER IN THE DERIVATIVE CQNLSE 

The discussion in section |TT] indicates that the loglog- 
normal divergence of the amplitude PDF is quite general 
for solitary waves of the CQNLS model. We now show 
that similar statistical behavior is exhibited by solitary 
waves of a second nonlinear wave model. We consider the 
derivative CQNLSE (DCQNLSE) with weak disorder in 
the linear gain/loss coefficient 
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where £(z) satisfies Eq. ([2]). In the context of non- 
linear optics e s is the self-steepening coefficient and 
ie s d t (\^\ 2 tp) describes the self-steepening effect. In the 
absence of the perturbation term ie£(z)i(), Eq. Jl 
possesses solitary wave solutions of the form tp a (t,z 
^ s (x) exp(ix), where 



{2-e s f3) 1 ' 2 p 



[(l-p 2 /^)V2 CO sh(2x) + l] 



1/2 ' 



(16) 



p = 17/(1 - e./3/2), p m = [4( £? - 3e2/16)/3]-V2, x = 
a + /3(t — y) + (r; 2 — /3 2 )z + g(x), and x — r)(t — y — 2(3z). 
In addition, the chirp g(x) is given by 



j(x) = p£sp m arctanh [B\ tanh(x)] , 



where the coefficient B± is 

^i-(i- P 2 A4) 1/2 



Bi = 



i + (i-p7j4)V2 



1/2 



(17) 



(18) 



These solitary wave solutions exist provided that e 9 > 
3e 2 /16 and p < p m . Notice that the solitary waves of 
the DCQNLSE are chirped and are thus fundamentally 
different from the solitary waves of the CQNLSE. 

In the presence of disorder in the linear gain both 77 
and randomly vary along the propagation. Thus, the 
dynamics is different from that observed in the CQNLSE 
case, where only 77 varies as a result of the disorder. 
Energy-balance considerations lead to an equation of the 
form 



dt\^\ 2 = 2e£(z) / dt\il>\ 



(19) 



In the first order adiabatic perturbation procedure we 
replace ip(t,z) with ip s (t,z) in Eq. (|19p and obtain 



d_ 

dz 



arctanh 



P 



= 2e£(z) arctanh ( — 



Denoting pd{z) = arctanh [p(z)/p r , 
satisfies the stochastic equation 



(20) 

we observe that pd 



dp d /dz = 2e£(z) Pd , (21) 

whose solution in Stratonovich's interpretation is: 

Pd(z)=p d (0)exp[2eX(z)]. (22) 

Therefore, the PDF of Pd{z) is lognormal. Changing vari- 
ables from p d (z) to p(z) we obtain that the PDF of p is 
given by 



exp {- In 2 [arctanh (p/p m ) /p d (0)} /{8De 2 z)} 
Z (8irDe 2 z) 1 / 2 p m (l-p 2 /p 2 m )&rctanh{p/p m ) ' 



(23) 



for < p < p m and F d {p) = elsewhere. In Appendix 
[B]we obtain the same expression for F d {p) by using the 



Fokker-Planck approach. The PDF Fd(p) has exactly 
the same form as the PDF of 77 in systems described by 
the CQNLSE with disorder in the linear gain/loss coeffi- 
cient. (Compare Eq. with Eq. (6) in Ref. In 
particular, Fd(p) exhibits loglognormal divergence in the 
vicinity of p rn : 



F d (p)\ p ^ pm c[(87rDe 2 z) 1 / 2 Sp\\n[Sp/(2p m )}\y 1 
x exp {- ln 2 [- In [Sp/(2p m )\ /(2p d (0))] /(SDe 2 z)}, (24) 



where Sp = p rn — p and < Sp/p m <C 1. In Appendix [Cl 
we calculate Fd(p) by applying Ito's interpretation [58j, 
[59j to the stochastic equation (f^T|) satisfied by p d . We 
show that in this case as well Fd(p) exhibits loglognormal 
divergence near p m . 

To validate our theoretical predictions we performed 
Monte Carlo simulations with Eq. (JT5J) with about 
7.8 x 10 4 disorder realizations. We used an initial condi- 
tion in the form of the solitary wave ip s with 77(0) = 1, 
/3(0) = 0, y(0) = 0, and a(0) = and considered the pa- 
rameter values D = 3, e q — 0.7, e s — 0.8, and e = 0.05. 
For these values, p(Q) = 1 and p m ~ 1.13715. The simu- 
lations were carried out up to a final distance Zf = 11, for 
which the disorder strength is 4De 2 z/ = 0.33. Equation 
(fT5|) was integrated by employing the split-step method 
with periodic boundary conditions and with the same nu- 
merical scheme as described in Sec. |TTJ The size of the 
computational domain was taken to be —100 < t < 100, 
and the t-step and z-step were taken as At = 0.01 and 
Az = 0.001, respectively. 

The PDF Fd(p) at z = 10 obtained in the numerical 
simulations is shown in Fig. Q] together with the theoret- 
ical prediction. The numerically obtained PDF clearly 
exhibits divergence in the vicinity of p m and the overall 
agreement between theory and simulation is good. To 
further check the behavior of the PDF in the vicinity 
of p m Fig. [5] shows a blowup of the data in the re- 



gion p 



< 



p m . Reasonable agreement between theory 
and simulations is observed. We attribute the differ- 
ences between the curves to the difficulties in obtain- 
ing an accurate measurement of the group velocity from 
the numerical data. As a further test for the asymp- 
totic behavior of Fd(p) near p m we employ the proce- 
dure for detecting loglognormal divergence that is out- 
lined in Appendix[XJ Following this procedure we present 
the graph of G d (Sp) versus gd(5p) in Fig. [6l where 
G d (6p) = {-H(87TDe 2 z) 1 / 2 8p\ln[5p/(2p m )}\F d (p)}} 1 / 2 
and g d (Sp) = (8De 2 z)- 1 / 2 ln{- ln[Sp/(2p m )}}. We ob- 
serve that the numerically obtained curve lies on a 
straight line with a slope 0.97, which is very close to 
the theoretically predicted value of 1. We therefore con- 
clude that the numerically obtained PDF of p does ex- 
hibit loglognormal divergence. 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
P 



2.0 2.4 2.8 3.2 

g d (5p) 



FIG. 4: The probability density function of p Fd(p) at z=10 
for D = 3, e q = 0.7, e s = 0.8, e = 0.05 and p(0) = 1. The 
squares represent the result of Monte Carlo simulations with 
Eq. (|15[) , while the solid curve corresponds to the analytic 
result given by Eq. (|23p . 



FIG. 6: Gd(Sp) vs gd(Sp) for the same parameters considered 
in Figs. [4] and [5] The squares represent the analytic result 
obtained with Eq. (|23[). while the circles stand for the nu- 
merical result. The solid and dashed lines are linear fits with 
slopes 1.00 and 0.97, respectively. 




1.02 1.04 1.06 1.08 1.10 1.12 1.14 
P 

FIG. 5: Blowup of the data shown in Fig. [4] in the vicin- 
ity of p m . The dashed line corresponds to the asymptotic 
loglognormal PDF given by Eq. (|24[1 , 



IV. DISCUSSION 

Here we discuss the underlying reason for the loglog- 
normal divergence of the 77-PDF due to disorder in the 
linear/nonlinear gain/loss. In addition, we briefly discuss 
a KdV-type of model where the loglognormal divergence 
of soliton parameters can potentially be observed. 



Using the adiabatic perturbation method we obtain the 
following equation for the dynamics of 77: 



d V 

-J- = 2er M 

az 



—5- I arctanh 



(-) 



(26) 



Even though Eq. (|26[) can be solved analytically, it is 
instructive to consider its asymptotic approximation for 
77 < r] m . Denoting Srj = r) m — r\ and expanding both sides 
of Eq. (f2"6"|) about r\ m while keeping terms up to 0(Sij) 
we obtain 



— ~ 2e<5r/ln ( — 
dz \ 2r) m 



Integrating Eq. (|27p over z we arrive at 



Sr](z) ~ 2r] m exp — C(0)e 



2ez 



(27) 



(28) 



where C(0) = - ln[^(0)/(2r? m )] > 0. We therefore ob- 
serve that in this case Srj decays to super-exponentially 
with increasing z. It is this super-exponential approach 
of 77 to rj m that leads to the loglognormal divergence of 
F(rj). Indeed, for the CQNLSE with disorder in the lin- 
ear gain/loss coefficient one obtains a similar equation 
for Sr](z) with z replaced by X{z) on the right hand side. 
As a result, X(z) is related to Srj(z) via 



A. Loglognormal divergence of F(rf) and 
super-exponential decay of <5?7 to 

Consider, for example, the CQNLSE with determinis- 
tic linear gain: 

id z *p + <9 t 2 V + 2|-0| 2 V> - e,MV = ^ (25) 



X(z) ~ In 



/(2e), (29) 



which describes loglognormal divergence of the ?7-PDF 
near r] m . A similar result holds for the CQNLSE with 
disorder in the cubic nonlinear gain/loss coefficient [see 
Eq. CGP]. 
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B. The extended Korteweg-de Vries equation 

A different type of nonlinear wave equation that pos- 
sesses fiat-top solitary wave solutions is the following; ex- 
tended Korteweg-de Vries (eKdV) equation 3, 6, 48]: 



dtu + 6u(l — u)d z u + d^u = 0. 



(30) 



Note that Eq. (J30J) is integrable [6|,|48|. In the context of 
interfacial waves in two-layer systems and stratified shear 
flow in the ocean u represents the vertical displacement, 
z is the horizontal coordinate, and t is time. The solitary 
wave solutions of Eq. (|30j) take the form [y] 



u s {z,t) 



An 2 



(1 - K 2 /^) 1 / 2 cosh(2a;) + l' 



(31) 



where x = k(z — An 2 t), n m — 1/2, and the parameter k 
characterizes the soliton amplitude and group velocity. 

Since perturbative linear gain and loss terms are quite 
common in KdV models of wave motion in the ocean 
[EH, [62j , it is interesting to study the situation where ran- 
domness is present. We therefore consider the following 
perturbed eKdV equation: 



where 



d t u + 6u(l - u)d z u + d^u = e£(t)u, (32) 



(£(<)> =0, mm)=DS(t-t'). (33) 



Employing mass-balance considerations we obtain that 
the dynamics of k is described by an equation similar to 
Eq. (|2"0)) . Based on this observation one would expect the 
k-PDF to exhibit loglognormal divergence in the vicinity 
of n m . It should be pointed out, however, that for models 
of the KdV type radiative effects are more significant 
[62| . and as a result, a perturbative calculation that takes 
these effects into account is required. We therefore defer 
the full analysis of this case to a future publication. 



weak multiplicative dissipative disorder. Furthermore, 
we showed that this divergence can be explained by the 
super-exponential approach of r\ to r\ m in the correspond- 
ing deterministic model with weak linear/nonlinear gain. 

Next we considered propagation of solitary waves in 
the presence of weak disorder in the linear gain/loss 
in systems described by the derivative cubic-quintic 
nonlinear Schrodinger equation (DCQNLSE). The soli- 
tary waves of the DCQNLSE are chirped and are thus 
fundamentally different from the solitary waves of the 
CQNLSE. As a result of the chirp, in the presence of dis- 
order in the linear gain/loss both the amplitude r\ and the 
group velocity (3 vary randomly along the propagation. 
We therefore studied the PDF of the parameter p, where 
p = 77/(1 — e s /3/2) and e s is the self-steepening coeffi- 
cient. We found that Fd(p) is loglognormally divergent 
when p is near its maximum value p m , i.e., the p-PDF of 
the corresponding flat-top solitary waves exhibits loglog- 
normal divergence. Moreover, we showed that the same 
divergence is observed in both Stratonovich's interpreta- 
tion and Ito's interpretation of the linear stochastic per- 
turbation term, thus illustrating another feature of the 
statistics that appears to be quite general. 



APPENDIX A: THE G VS g METHOD 

Here we give the details behind the G vs g method 
that is used to analyze the divergence of the PDFs of r\ 
and p near their maximum possible values. As a spe- 
cific example we consider the case of the CQNLSE in the 
presence of disorder in the cubic gain/loss coefficient. In 
Sec. HI we obtained the following approximate analytic 
expression for F c (rf) near r/ m : 



F c(v)\ v < v „ 



In 



e 2 Sr] 
2r)m 



xexp 



In 2 [-In [(e 2 5 V )/(2 Vm )] /g] 



(Al) 



V. CONCLUSIONS 

We studied the statistics of flat-top solitary wave pa- 
rameters in the presence of weak dissipative multiplica- 
tive disorder. We started by considering propagation of 
solitary waves of the cubic-quintic nonlinear Schrodinger 
equation (CQNLSE) in the presence of disorder in the 
cubic nonlinear gain/loss. We found that the ampli- 
tude PDF exhibits loglognormal divergence in the vicin- 
ity of the maximum possible amplitude ry m . Since soli- 
tary waves with rj values near r\ m have a typical table-top 
shape we conclude that the amplitude PDF of flat-top 
solitary waves is loglognormally divergent. This find- 
ing combined with similar findings in Ref. [54| for the 
case of disorder in the linear gain/loss coefficient indi- 
cates that loglognormal divergence is quite ubiquitous for 
flat-top solitary waves of the CQNLSE in the presence of 



The problem we want to address is how to verify that 
the PDF obtained in the simulations satisfies the loglog- 
normal divergence described by Eq. (jAljl . In particular, 
we need a method that will allow us to ignore the co- 
efficient c, which cannot be found from the numerical 
data. Furthermore, we would like to find a mapping that 
"stretches" the small ^-neighborhood of r/ m into a wider 
interval. To address these issues we rewrite Eq. (|Alj) in 
the form 



In 



(32^ De 2 zf/ 2 V 2 m 5 V 



In 



/ e 2 8-q 
\2?7 m 



ln 2 {-ln[(e 2 ^)/(2r ?m )] /g} 
32De 2 V iz 



(A2) 



where it is understood that F c (rj) is calculated near r] m . 
Multiplying by -1 and taking the square root we arrive 
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at 



-In 



(32wDe 2 z)^ 2 V 2 m Sr] 



In 



e 2 Sr] 
2r) m 



ln{-ln [(e 2 5 V )/(2 Vm )] } - ln(c) 

(32i^*) 1/a 
We now define G c (5r]) and g c (Sv) as 

Gc(^) = 



In 



(32nDe 2 z) 1 / 2 7 1 2 m Sr 1 



In 



and 



9c(Srj) 



ln{-ln[( e 2 ^)/(277 m )]} 
(32De^ V f n zf 2 



1/2 



(A3) 



1/2 



(A4) 



(A5) 



Using these definitions we observe that when the statis- 
tics of r\ is described by Eq. (|A1[) the graph of G vs g is 
a straight line with a slope 1, independent of the value 
of c. A similar conclusion (with slightly different expres- 
sions for G and g) holds in the case of DCQNLSE with 
disorder in the linear gain/loss coefficient. 



for < p < p m and F d (p) — elsewhere. Equation 
(|B5|) is the same as Eq. (I2U1) in section IIIII A similar 
calculation based on the Fokker-Planck approach leads 
to Eq. Uni) for the r]-PDF for the CQNLSE with disorder 
in the cubic gain/loss coefficient. 



APPENDIX C: CALCULATION OF F d {p) IN ITO'S 
INTERPRETATION 



Consider the DCQNLSE with disorder in the linear 
gain/loss coefficient. We now obtain the PDF of p by 
employing ito's interpretation to equation (|20[) , and show 
that this PDF exhibits loglognormal divergence in the 
vicinity of p m . The solution of the equivalent equation 
(l2"Tj) in Ito's interpretation is [HJ 

Pd{z) = p d (0) exp [2eX(z) - 2e 2 z] , (CI) 

where the PDF of X(z) is given by Eq. ©. Therefore, 
X(z) is related to p(z) via 



APPENDIX B: FOKKER-PLANCK APPROACH X(z) 
FOR CALCULATION OF THE PDFS 



hi 



arctanh 



p( z ) 

Pm 



/MO) 



2e z z 



(C2) 



In this Appendix we demonstrate that the expressions 
for the PDF of the solitary wave parameters that were ob- 
tained in sectionsHTland lllll bv solving the Langevin equa- 
tion can also be obtained within the framework of the 
Fokker-Planck approach. As a specific example we con- 
sider the DCQNLSE with disorder in the linear gain/loss 
coefficient and work with Stratonovich's interpretation of 
Eq. (|2ip . The corresponding Fokker-Planck equation for 
the PDF of p d , H{p d , z), is li, [H: 



d z H = 2De 2 d Pd [ Pd d Pd (p d H)) 
Changing variable to w = In p d we arrive at 
d z H = 2De 2 dlH, 

where 

H(w, z) = p d {w)H{p d (w), z). 



(B2) 



The solution of Eq. (|B2j) with the initial condition 
H (w, z) = S(w — w(0)) is 



H(w, z) 



Changing variable from 
and (IB4I) we obtain 



exp {-[w - w{0)] 2 /(8De 2 z)] 
(87r J De 2 z) 1 /2 ■ 

to p while using Eqs. 



(B4) 



exp {- ln 2 [arctanh (p/ Pm ) /Pd(0)} /(8De 2 z)} 
(%nDe 2 z) l / 2 p m (l~p 2 /p 2 Il ) arctanh (p/p m ) 



(B5) 



Changing variables from X(z) to p(z) we obtain the PDF 
of p in Ito's interpretation: 



F^ I \p)^\{87:De 2 z) 1 ^p m (l-p 2 /p 2 m ) arctanh (p/ Pm ) 



x exp 



{- [ln[arctanh (p/p m ) /p d (Q)} + 2e 2 z] 2 /(8De 2 z)} . 

(C3) 



(Bl) In the vicinity of p m Eq. (|C2[) can be approximated by 



X(z) 



hi 



2 Pd (0) 



■In 



5p{z) 



2p n 



+ 2e 2 z )■ . (C4) 



(B3) Consequently, (p) is given by 



F, 
exp 



f\p) ^U%nDe 2 z) 1 ' 2 8p\\n[8p/{2p m )]\Y\ 
{- [ln[- ln[8p/(2p m )]/(2 Pd (0))] + 2e 2 z} 2 /(8De 2 z)} 



(C5) 



which exhibits loglognormal divergence as p approaches 
p m . We therefore conclude that for the DCQNLSE with 
disorder in the linear gain coefficient both Stratonovich's 
interpretation and Ito's interpretation lead to loglognor- 
mal divergence of the PDF of p. 
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